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Can the third dredge-up extinguish hot-bottom burning in 

massive AGB stars? 
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Abstract. Marigo (2002) has highlighted the crucial importance of molecular opacities in modelling the evolution of AGB 
stars at varying surface C/O ratio. In particular, it has been shown the large inadequacy of solar-scaled opacities when applied 
to models of carbon stars, and hence the need for correctly coupling the molecular opacities to the current surface chemical 
composition of AGB stars. The aim of the present follow-up study is to investigate the effects of variable molecular opacities 
on the evolutionary properties of luminous AGB stars with massive envelopes, i.e. with initial masses from » 3.5M Q up to 
5-8 M Q , which are predicted to experience both the third dredge-up and hot-bottom burning. It is found that if the dredge-up of 
carbon is efficient enough to lead to an early transition from C/O < 1 to C/O > 1, then hot-bottom burning may be weakened, 
extinguished, or even prevented. The physical conditions for this occurrence are analysed and a few theoretical and observational 
implications are discussed. Importantly, it is found that the inclusion of variable molecular opacities could significantly change 
the current predictions for the chemical yields contributed by intermediate-mass AGB stars, with M = 3.5 - 4.0 M Q , that make 
as much as ~ 30 - 50 % of all stars expected to undergo hot-bottom burning. 



1. Introduction 

During the Thermally-Pulsating Asymptotic Giant Branch 
(TP-AGB) phase, intermediate-mass stars (with initial masses 
3.5M $ M < M ap = 5 - 8M , depending on metallicity and 
model details) experience a particularly rich nucleosynthesis, 
whose products are exposed to the surface via two main pro- 
cesses, namely: the third dredge-up, and hot-bottom burning 
(hereinafter also HBB). 

The former injects into the envelope elements like helium 
and other species produced during thermal pulses, essentially 
primary carbon and s-process elements (Iben & Renzini 1983; 
Busso et al. 1999). The formation of carbon (C-) stars of N type 
(with surface C/0> 1 and not belonging to binary systems) is 
naturally explained in this way. 

The latter process corresponds to nuclear H-burning that 
extends from the radiative shell into the deepest and hottest 
layers of the convective envelope during the quiescent inter- 
pulse periods. The main nuclear reactions at work are those of 
the CNO-cycle, so that the surface abundances of helium and 
nitrogen are expected to increase. Heavier nuclei - like some 
isotopes of neon, sodium, aluminum, and magnesium - can be 
produced from the Ne-Na and Mg-Al cycles, whereas some 
amount of fresh lithium can also be synthesised and brought 
to the surface via the so-called beryllium-transport mechanism 
(Cameron & Fowler 1971), which requires suitable conditions, 
involving the temperature at the base of the convective enve- 
lope, and the time scales of nuclei and convective eddies (see 
Lattanzio & Wood 2003 for an overview of the HBB nucle- 
osynthesis). 



Actually, HBB has been invoked to explain (or to con- 
tribute explaining) a number of observational evidences, like 
for instance: the (almost) lack of visible luminous C-stars 
(Mboi < -6) and the concomitant deficiency of luminous M- 
stars in the Magellanic Clouds (e.g. Costa & Frogel 1996); 
the existence among these stars of objects (-6 > Mb i > -7) 
with marked enhancement of surface lithium abundance (e.g. 
Smith et al. 1995); the overabundance of nitrogen and helium 
typical of type I Planetary Nebulae and some AGB stars (e.g. 
Peimbert & Torres-Peimbert 1983; Smith & Lambert 1990); 
the low value of the 12 C/ 13 C ratio, close to the nuclear equi- 
librium value, characterizing the group of luminous J-type car- 
bon stars (e.g. Lambert et al. 1986); the Na-O and Mg-Al anti- 
correlations characterising the stars of galactic globular clus- 
ters in the framework of the primordial scenario (Gratton et 
al. 2001); the high isotopic ratios 26 A1/ 27 A1 displayed by some 
meteoritic oxide grains of pre-solar origin (Mowlavi & Meynet 
2000). 



Several evolutionary models of massive AGB stars with 
HBB (say with M ^ 3.5 - 4M , depending on metallic- 
ity) have been calculated in recent years (e.g. Boothroyd & 
Sackmann 1992; Blocker 1995; Forestini & Charbonnel 1997; 
Marigo 1998; Frost et al. 1998; Ventura et al. 1999; Karakas & 
Lattanzio 2003; Herwig 2004; Ventura & D'Antona. 2005ab). 
Compared to earlier studies, these calculations mark a net 
progress in respect to many aspects, like for instance the adop- 
tion of larger nuclear networks including many elemental and 
isotopic species, a more detailed description of convective tur- 



2 



P. Marigo: Extinction of Hot Bottom Burning by the Third Dredge-up 



bulence and mixing, and extensive calculations over wider 
ranges of initial stellar masses and metallicities. 

However, one standard choice common to all these mod- 
els is the use of tables of molecular opacities - for tempera- 
tures lower than ~ 10 000 K - that are strictly valid for solar- 
scaled abundances of elemental species heavier than helium 
(e.g. Alexander & Ferguson 1994). Then, for a fixed initial 
metallicity, interpolations are usually performed as a function 
of temperature, density, and hydrogen content. This implies 
that any change in the elemental and molecular concentrations 
(following the convective dredge-up and HBB) does not trans- 
late, as instead it should, into a change of the molecular opaci- 
ties adopted to calculate the AGB model atmospheres. 

This point of inconsistency has been clearly pointed out 
by Marigo (2002; see also Marigo et al. 2003b), who has con- 
structed a routine to estimate the molecular opacities for what- 
ever chemical composition of AGB envelopes. Typically, the 
photospheric opacities are dominated by an opacity bump due 
to H2O for C/O < 1, while as soon as C/O > 1 this peak disap- 
pears and a prominent opacity bump due to the CN molecules 
develops. 

The routine has been then incorporated into a static enve- 
lope model, which is part of a code developed to calculate the 
synthetic TP-AGB evolution (see Marigo et al. 1999 and refer- 
ences therein). By comparing the results of TP-AGB evolution- 
ary calculations obtained with fixed solar-scaled to those with 
variable molecular opacities it turns out that the surface C/O ra- 
tio plays a key-role in determining the evolutionary properties 
of AGB stars. 

In Marigo (2002) most attention was devoted to study the 
impact of variable opacities on the evolutionary behaviour of 
low-mass stars (with masses up to 2.5M Q ). It turns out that 
one of the major effects is the sudden and marked cooling of 
the evolutionary tracks in the H-R diagram as soon as C/O 
overcomes unity, that is at the transition from the oxygen-rich 
configuration (C/O < 1) of M-type stars to that of carbon-rich 
stars (C/O > 1). The attainment of lower effective temperatures 
causes, in turn, other important evolutionary consequences, 
namely: earlier onset of the super-wind; shorter duration of the 
TP-AGB and C-star phases, redder near-infrared colours, re- 
duced degree of carbon enrichment hence lower chemical yield 
and C/O ratios in carbon stars. Most of these predictions al- 
lows a better agreement with the observations of carbon stars. 
Moreover, the relatively low C/O ratios ($1) exhibited by disk 
planetary nebulae can be more easily explained with the aid 
of TP-AGB models including the third dredge-up and variable 
molecular opacities, as shown by Marigo et al. (2003a). 

Following the introduction of the new synthetic TP-AGB 
models into a population synthesis code Marigo et al. (2003b) 
have shown that the variable molecular opacities are the key 
to explain the prominent red tail of visible carbon stars seen 
in some near-infrared colour-magnitude (e.g. K vs. J - K) di- 
agrams, like those provided by the DENIS (Cioni et al. 2000) 
and 2MASS surveys (Nikolaev & Weinberg 2000) towards the 
Magellanic Clouds. 

The present paper pushes forward the analysis started with 
Marigo (2002) and Marigo et al. (2003b), by investigating the 
possible effects driven by variations of molecular opacities in 



the most massive AGB stars (with M <: 3.5M Q ), that experi- 
ence both the third dredge-up and hot-bottom burning. In par- 
ticular it will be shown that, if during the early stages of the 
TP-AGB phase the third dredge-up is able to produce C/0> 1 
in the envelope, then this may prevent, or weaken and eventu- 
ally extinguish the nuclear rates associated with HBB. 

It follows that the third dredge-up and HBB produce re- 
ciprocal interference. In fact, on one side, it has been long 
predicted (e.g. Renzini & Voli 1981; Boothroyd & Sackmann 
1992) that HBB can delay or even prevent the formation of 
high-luminosity C-stars, by destroying the dredged-up carbon 
in favour of nitrogen. On the other side, the present study 
shows, for the first time, that the third dredge-up itself may con- 
trast the development of HBB. This means that the two com- 
peting processes are intimately related one to each other, in a 
more complex manner than considered so far. 

The structure of the paper is as follows. Section [2] details 
the envelope model adopted in the present study. The sensitive- 
ness of luminous and massive stellar envelopes to increasing 
surface C/O ratio is first explored in Sect. [3]with the aid of en- 
velope integrations, and then better investigated in Sects. |4]and 
where a few synthetic TP-AGB calculations are presented 
for selected values of the stellar initial mass and metallicity. 
Next, in Sect. [6] the analysis is focused on the impact of vari- 
able molecular opacities in changing the predicted evolutionary 
properties of a particular group of stars, i.e. the lithium-rich 
carbon stars, while Sect. [7] weighs the significance of the new 
results in the context of the theory of stellar populations. The 
reader is made aware of a few cautionary remarks on possible 
limits of the present work in Sect. [8] whereas a closing sum- 
mary is provided in Sect. [9] 

2. Envelope integrations 

The calculations presented in this paper are based on numerical 
integrations of a static envelope model, which is also a key- 
ingredient of the synthetic TP-AGB model already developed 
by Marigo et al. (1996, 1998, 1999). 

Let us just briefly outline here the basic scheme. Envelope 
integrations are carried out from the photosphere down to 
mass-coordinate of the H-exhausted core. The determination 
of the unknown functions r, P r , T r , L r across the envelope is 
obtained by solving the stellar structure equations, once four 
boundary conditions are specified. Two conditions naturally 
derive from the integration of the photospheric equations for 
T and P down to the bottom of the photosphere (Kippenhahn 
et al. 1967). The remaining two conditions involve the radius 
and the local energy flux at the mass-coordinate of the He-core 
(see Marigo et al. (1998), Marigo (1998) for all the details). 

Given the mass of the core, and the mass and chemical com- 
position of the envelope, this procedure allows to single out the 
main envelope parameters, namely the luminosity L and the ef- 
fective temperature T e ff . 

It is important to remark that our model accounts for nu- 
clear energy sources possibly operating in the deepest en- 
velope layers. The common prescription of envelope models 
L r = constant = L is therefore replaced with the equation of 
energy balance dL r /dM r = e, (with usual meaning of the quan- 
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Fig. 1. From top to bottom: predicted behaviour of the temper- 
ature at the base of the convective envelope, effective tempera- 
ture, and luminosity at varying surface C/O ratio. Carbon abun- 
dance (by number) is made to increase in steps of 0.05, while 
that of oxygen is kept constant. The curves are the results of 
envelope integrations, assuming a total stellar mass of 4.5 M , 
a core mass of 0.98 M , an initial metallicity Z = 0.004,, and 
adopting either k^ x or k vm . See text for more explanation 



tities). In this way we can consistently follow hot-bottom burn- 
ing in the more massive AGB stars, with respect to both the 
CNO-cycle nucleosynthesis and the related over-luminosity ef- 
fect that makes these models more luminous than predicted by 
the core mass-luminosity relation for the same M c (see Marigo 
et al. 1998, Marigo 1998 for a detailed discussion). 

The basic physical inputs employed in our envelope model 
are as follows: Nuclear reaction rates for the p-p chains and 
CNO-bicycle are taken from the compilation by Caughlan & 
Fowler (1988), the screening factors are given by Graboske et 
al. (1973). The adopted mixing-length parameter is a = 1.68. 

As far as the gas opacities are concerned we use those cal- 
culated by Iglesias & Rogers (1996) (OPAL) at high tempera- 
tures (T > 10 4 K). At lower temperatures (T < 10 4 K) we adopt 
the opacity routine developed by Marigo (2002) choosing be- 
tween two alternatives for the adopted chemical composition: 

- a solar-scaled mixture of elemental species at any metal- 
licity, as in Alexander & Ferguson (1994) opacity tables 
(hereinafter also ksx); 

- the current mixture of chemical elements as it results from 
the nucleosynthesis and mixing events occurring during the 
TP- AGB evolution of the model under consideration (here- 
inafter also K V ar)- 



3. Sensitiveness of extended AGB envelopes to 
the C/O ratio 

We will now investigate the effects produced by an efficient 
dredge-up of carbon on the envelope structure of massive and 
luminous AGB models, that are also predicted to experience 
HBB. 

To this aim, Fig. Q] shows the results of some test calcula- 
tions, carried out by means of envelope integrations with either 
/Cfi x or K var (see Sect.EJ. We assume a core mass M c = 0.98 M Q , 
a total stellar mass M = 4.5 M©, and an initial chemical com- 
position of the envelope [Zj = 0.004, F; = 0.24], with Z, Y l7 
X[ = 1 - Y{ — Zj being the abundances (in mass fraction) of 
metals, helium and hydrogen, respectively. 

Initially, the elemental abundances of all species heavier 
than helium are assumed to be solar-scaled according to com- 
pilation of solar abundances by Grevesse & Noels (1993). It 
follows that the starting C/O ratio in the envelope is ~ 0.48. 
Then, envelope integrations are performed at varying C/O ratio, 
which is made increase from below to above unity by adding 
carbon ( 12 C) while keeping oxygen fixed to its solar-scaled 
abundance. We specify that in all models discussed in the paper 
the C/O ratio is calculated 

C _ X( 12 C)/12 + X( 13 C)/13 

O ~ X( 16 0)/16 + X( 17 0)/17 + X( ls O)/18 ' " ' 

where the Xs in the right-side member denote the elemen- 
tal abundances of the C and O isotopes (mass fractions). The 
excess of carbon over its solar-scaled value, AC = X( 12 C) - 
X( 12 C) B , is compensated by reducing the hydrogen abundance 
by the same amount, that is X = X\ - AC, which also implies 
that Y = Yi. 

FigureQ]displays the behaviour of relevant quantities - i.e. 
temperature at the base of the convective envelope 7b ce , effec- 
tive temperature r e ff, and stellar luminosity L - as a function of 
C/O. 

The first point to be noticed is that, while in the case 
the variables show a smooth and regular trend over the whole 
C/O range, a marked slope change shows up at the transition 
from C/O < 1 to C/O > 1 in the /c var curves. Second, with the 
/(•fix assumption the variations of 7bce, TW, and L are globally 
smaller, and they may even proceed to the opposite direction 
compared to the /c V ar case, as is for the the surface luminosity. 
These differences can be explained as follows. 

On one hand, the steady brightening of the model with 
is essentially due to a cycling effect, that is more and more 
carbon is progressively made available to enter the CN cycle 
at the base of the convective envelope, thus increasing the ef- 
ficiency of nuclear reactions. Despite the modest decrease of 
Tbce, which would tend to lower the nuclear rates, this carbon 
cycling effect dominates in determining the behaviour of L. In 
the outermost atmospheric layers the effect of the increased 
carbon abundance on the temperature stratification becomes al- 
most negligible so that remains practically constant. 

On the other hand, the inflection points at C/0~ 1, that 
are present in all the three curves with /c var , are caused by the 
sudden change of the dominant opacity sources at low temper- 
atures (T < 6000 K), reflecting, in turn, the abrupt alteration 
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of the molecular concentrations (see Marigo 2002). In fact, the 
development of the large CN opacity bump for C/O > 1 causes 
a marked decrease of both T e ff and 7b ce - In this case, the cool- 
ing effect of the envelope base prevails over the injection of 
new carbon into the CN cycle, so that the efficiency of nuclear 
energy generation (ecNO K r 15 ~ 20 ) becomes lower with conse- 
quent decrease of the surface luminosity. 

In summary, from these simple test calculations we al- 
ready derive the indication that molecular opacities - coupled 
to the actual chemical abundances in the envelope - could sig- 
nificantly influence the structure and nucleosynthesis of mas- 
sive AGB stars by affecting the efficiency of HBB. Anyhow, it 
should be noticed that these preliminary results indicate how 
the structure of a stationary envelope reacts to changes in its 
chemical abundances, hence in molecular opacities, while all 
other parameters (e.g. core mass, total stellar mass) are kept 
fixed. The next necessary step is to analyse the impact on the 
evolutionary properties of massive AGB stars, which is done in 
the next section. 

4. TP-AGB evolutionary calculations 

We have calculated the TP-AGB evolution of a few stel- 
lar models experiencing both the third dredge-up and HBB. 
Calculations are performed with the aid of the synthetic TP- 
AGB model as recently revised by Marigo & Girardi (2006, 
and references therein), to whom we refer for all details (see 
also Sect. |2j. 

Suffice it here to recall that the initial conditions at the first 
thermal pulse - i.e. total mass, core mass, luminosity, surface 
chemical composition - are taken from Girardi et al. (2000). 
The third dredge-up is parametrized in terms of three quanti- 
ties, namely: the minimum core mass required for the occur- 
rence of convective dredge-up, M™ n ; the dredge-up efficiency, 
A; and the chemical composition of the dredged-up material. 

In practice, following the indications from full AGB calcu- 
lations (Karakas et al. 2002), in all models here considered the 
third dredge-up is assumed i) to occur since the first thermal 
pulse, and ii) to be quite efficient, typically with A « 0.9 - 1 .0. 
The adopted inter-shell abundances depend on the core mass 
growth as described in Marigo & Girardi (2006), reaching 
standard maximum values [X(He) 0.73, X( 12 C) - 0.25, 
X( l6 0) 0.02]. 

As already mentioned in Sect. [2] HBB is followed in de- 
tail with the aid of a complete envelope model, including the 
most important reactions rates of the CNO bi-cycle. The over- 
luminosity effect due to HBB and the consequent break-down 
of the core mass-luminosity relation are taken into account. 
Mass loss is described with the formalism by Vassiliadis & 
Wood (1993), under the assumption that variable AGB stars 
are fundamental-mode pulsators and applying the delay in the 
onset of the super-wind for M ^ 2.5 M (their equation 5). 

Let us first consider the case of a (3.5 M Q ,Z = 0.004) 
model. Figure |2]displays the evolution of characteristic quanti- 
ties - i.e. luminosity, effective temperature, temperature at the 
base of the convective envelope, and a few surface elemental 
ratios - from the first thermal pulse up to the termination of the 
TP-AGB phase, marked by the complete ejection of the enve- 



lope. The results are shown for both opacity prescriptions, Kg x 
(left panel) and A: var (right panel). 

From the comparison between the two models we notice 
that significant differences show up as soon as the surface C/O 
ratio overcomes unity as a consequence of the third dredge-up. 

First of all, we recover the same photospheric cooling ef- 
fect found for lower mass AGB models (say, with M < 3.5M ) 
that undergo the transition to the C-star domain, as fully dis- 
cussed by Marigo (2002) and Marigo et al. (2003). In fact, a 
sudden jump to lower r e ff is seen to occur in the /c val model at 
C/0<; 1, while a smooth temporal behaviour of T e ff still char- 
acterises the evolution with k^. The attainment of lower by 
the /c va , model yields longer pulsation periods, hence favouring 
an earlier onset of the super-wind regime of mass loss. As a 
consequence, the duration of the TP-AGB phase turns out to be 
reduced by almost a factor of 2. 

In addition to this, the introduction of variable molecular 
opacities is expected to impact significantly also on the nucle- 
osynthesis associated with HBB. This is the combined result of 
enhanced mass-loss rates and somewhat lower temperatures at 
the base of the convective envelope (see Sect. [3). Both factors 
concur to weaken the efficiency of nuclear reactions. 

As for the evolution of the C/O ratio, the following aspects 
should be noticed. As we see in Fig. [2] the model first be- 
comes carbon-rich (C/O > 1) because of the third dredge-up, 
and later it is converted back as oxygen-rich (C/O < 1) by the 
CN-cycle operating in the inner layers of the convective en- 
velope. This configuration is held up to the end of the AGB 
phase. On the contrary, in the /c val case, once the surface carbon 
abundance overcomes that of oxygen, the C/O ratio keeps on 
increasing so that at the termination of the TP-AGB phase the 
model is still carbon-rich. 

The fact that in the K val model HBB is effectively extin- 
guished soon after the transition to the C-star domain is clearly 
indicated by the behaviour of other surface abundances, which 
are displayed in the bottom panels ofFig.|2] The 12 C/ 13 C ratio, 
for instance, is prevented from attaining nuclear equilibrium 
(— 3 — 4), which is instead reached and maintained for long by 
the /ffi x model. More remarkably, the N/O ratio remains prac- 
tically constant over the entire TP-AGB evolution with *c var , 
implying that the CN cycle has never reached the necessary 
thermodynamical conditions (i.e. temperature and density) to 
convert carbon into nitrogen. Some nitrogen production is, in- 
stead, predicted for the model with k^ x , as shown by the rising 
part of the N/O curve. 

In summary, we expect that the efficiency of HBB may 
be reduced even significantly if, during the early stages of its 
TP-AGB evolution, a massive AGB star experiences efficient 
carbon dredge-up such as to become a carbon star. Then, the 
enhanced molecular opacities, associated with a carbon-rich 
chemical composition, importantly affect the temperature strat- 
ification of the layers above the core. A cooling effect takes 
place both at base of the convective envelope, i.e. lowering 7b Ce , 
and at the atmosphere, i.e. lowering T e ^. The way the two fac- 
tors concur to weaken, and possibly to extinguish HBB, is dif- 
ferent. The decrease of 7b ce has, clearly, a direct effect on the 
nuclear reaction rates of HBB, while the decrease of T e $ affects 
the HBB efficiency by altering the evolutionary path of the TP- 
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Fig. 2. From top to bottom: temporal evolution of luminosity, effective temperature, temperature at the base of the convective 
envelope, surface C/O, 12 C/ 13 C, and N/O ratios during the TP-AGB phase of a (M; = 3.5, Z = 0.004). Calculations are carried 
out from the first thermal pulse up to the ejection of the whole envelope according to Vassiliadis & Wood (1993) formalism for 
mass loss 



AGB model. In fact, a lower effective temperature favours the 
attainment of larger mass-loss rates, hence an earlier extinction 
of HBB due to the reduction of the envelope mass by stellar 
winds. 

It is clear that the prevailing effect between the two depends 
on several parameters such as: mass loss and its dependence on 
J e ff, dredge-up efficiency, stellar initial mass and metallicity. In 
Sect.|5]we will discuss exemplifying results for various choices 
of the stellar mass and metallicity. 



5. Dependence on mass and metallicity 

The interference of the third dredge-up against the efficiency of 
HBB develops to a different extent depending on several fac- 
tors, such as initial stellar mass, dredge-up characteristics (ef- 
ficiency A and chemical composition of the inter-shell), mass 
loss prescription, and metallicity. 

In this study we limit to explore the sensitiveness of the 
results on initial stellar mass and metallicity. A few cases are 
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M. = 4.0 Z = 0.008 
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Fig. 3. Predicted temporal evolution of surface elemental ratios (by number) during the whole TP-AGB evolution of a few 
selected models with varying initial stellar mass and metallicity, as indicated. Results are displayed for both choices of molecular 
opacities, Arg x and /c var . Both the third dredge-up and HBB are at work in these models 



presented in Fig. [3] that shows the time evolution of the surface 
chemical ratios C/O, 12 C/ 13 C, and N/O, adopting both k^ x and 
Kvar prescriptions. 

The three ratios give information on the strength of HBB 
and its competition with the third dredge-up, the behaviour of 
C/O showing directly if a star should appear either as an M- 
type or as a C-type giant, where the trends of 12 C/ 13 C and N/O 
measure the efficiency of HBB to activate the first nuclear re- 
actions of the CN cycle (via the conversion of 12 C to 13 C), and 



subsequently to produce primary N at the expense of the newly 
dredged-up C. 

The choice of the initial masses is determined by the re- 
quirement that, in the /cg x case, the models actually make an 
early transition to the C-star domain, which is the necessary 
condition to investigate a possible opacity effect in the cor- 
responding K yM calculations. In practice, this means to con- 
sider models experiencing a moderate HBB, that is with initial 
masses in the lowest part of the typical mass interval for its oc- 
currence, i.e. Mi ~ 3.5M Q - 4.0M o , depending on metallicity. 
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Fig. 4. Temperature, 7b ce , at the base of the convective envelope 
as a function of time during the TP-AGB evolution of stellar 
models with initial mass and metallicity as indicated. For the 
sake of simplicity only the quiescent pre-fiash stages are con- 
sidered in the construction of the plot. Results are shown for 
both /Cfi x (left panels) and K wm - (right panels) prescriptions. Thick 
lines correspond to C/O > 1 at the surface. The hatched strip 
is defined by the temperature conditions required to produce 
lithium-rich carbon stars according to Ventura et al. (1999). See 
text for more explanation 

From the inspection of Fig. [3] it turns out that the discus- 
sion made in Sect. 0] in respect to the (Mj = 3.5, Z = 0.004) 
model applies also to the (Mj = 4.0, Z = 0.019) and (Mj = 
4.0, Z = 0.008) models: The adoption of variable molecular 
opacities leads to an effective extinction of HBB shortly after 
the attainment of C/O > 1 . 

Different is the case of the (Mj = 3.5, Z = 0.001) model 
(bottom panel), with the lowest metallicity here considered. 
The differences of the results between k^ x and /c val are just mi- 
nor and both models experience efficient HBB despite the fact 
that they become carbon-rich during the initial stages of their 
TP-AGB evolution. This can be understood by considering that 
at decreasing metallicity molecular concentrations get lower 
and lower because of i) the lower abundances of the involved 
atoms, and ii) the strong temperature dependence of molecu- 
lar formation on temperature. In extremely metal-poor atmo- 
spheres the temperatures may be so high that the molecules 
cannot be even assembled. 

6. The evolution of lithium-rich carbon stars 

As shown in the previous sections, the opacity effect driven by 
the third dredge-up is expected to be important for AGB stars 
that both are carbon rich and experience HBB. Interestingly, 
there is one group of observed objects that could meet these 



conditions, namely the luminous lithium-rich carbon stars. In 
fact, among the most luminous carbon stars detected in the 
Magellanic Clouds a few possess a significant enhancement of 
lithium abundance in the magnitude range -5.5 ^ M\, \ S -6.5 
(e.g. Smith & Lambert 1989, 1990). 

These stars are currently explained considering the com- 
bined effect of the the third dredge-up for the surface enrich- 
ment of carbon, and the Cameron & Fowler (1971) mecha- 
nism, operating under HBB conditions, for the synthesis and 
temporary survival of lithium at the photosphere. Ventura et al. 
(1999) address this question with the aid of full AGB evolution- 
ary models (with initial masses 3.3< M/M Q < 4.0). They show 
that the chemical peculiarities of the high-luminosity lithium- 
rich C-stars can be reproduced only under strict structural con- 
ditions for the occurrence of HBB. Specifically, the tempera- 
ture at the base of the convective envelope should vary within a 
narrow range, that is 4 x 10 7 ^ 7b ce ^ 6.5 X 10 7 . The lower limit 
is the minimum temperature required to activate the Cameron 
& Fowler mechanism for lithium production, the upper limit 
corresponds to the maximum temperature to preserve carbon 
against an efficient destruction by the CNO-cycle. This temper- 
ature requirement would be fulfilled by AGB stars with initial 
masses and luminosities within small ranges. Just in virtue of 
the predicted narrow luminosity domain these stars might be 
also used as independent distance-indicators, as suggested by 
Ventura et al. (1999). 

In the light of the analysis just developed we can already 
envisage that the introduction of variable molecular opacities 
(not included in Ventura et al. (1999) calculations) may al- 
ter some of those findings. It is not among the scopes of the 
present paper to perform a detailed analysis of all these issues, 
however our test calculations show that the impact of variable 
molecular opacities may be important for models of lithium- 
and carbon-rich stars and future detailed models should include 
this physics. 

Figure H] helps to illustrate this point. The evolution of T\, ce 
is shown as a function of time during the TP-AGB phase of the 
same intermediate-mass models already discussed in the previ- 
ous sections, assuming either or k vm . Following Ventura et 
al. (1999) the strict temperature conditions for the formation of 
lithium-rich carbon stars define the hatched rectangular region. 

A systematic difference affects the results obtained with 
variable molecular opacities. While in models with /Cfj x the C- 
rich and Li-rich phase is characterised by a steady increase of 
Tbce which keeps on for some time even after the subsequent 
conversion to an oxygen-rich configuration, in models with /c var 
the same phase always intercept the maximum of r^e, after 
which an abrupt drop is expected as a consequence of the onset 
of the super-wind. 

As a consequence, we expect significant changes in the 
predicted properties of these stars, such as the duration, the 
luminosity range, and the maximum Li production of this 
chemically-peculiar phase. Important differences should result 
also in respect to the chemical yields from these stars. In fact 
the latter group of models with /c var is expected to enrich the 
ISM with newly synthesised carbon and lithium, while the for- 
mer models with /Cfj x should eject primary nitrogen and lithium. 



8 



P. Marigo: Extinction of Hot Bottom Burning by the Third Dredge-up 



A quantitative prediction of the yields from these stars cannot 
be given at this stage and it is postponed to a future work. 

7. Weighing the impact of the new results 

As already discuss in Sect. [5] the introduction of variable 
molecular opacities may importantly affect AGB models expe- 
riencing both the third dredge-up and HBB, within a relatively 
narrow range of initial stellar masses, roughly from 3.0 up to 
4.0 M Q , depending on metallicity and model assumptions. 

To assess the significance of the new results we will exploit 
a few basic relations from the theory of stellar populations (e.g. 
Tinsley 1980). In a composite stellar population like a galaxy, 
the expected number of TP-AGB stars with initial masses in 
the interval [m, m + dm] is expressed by 

N{m) dm oc <p(m) iJ/[Tg — t(mj] fTP-AGB('«) dm (2) 

where </>(w) is the initial mass function (IMF), i//[Tc - t(m)] is 
the star formation rate (SFR) evaluated at the time of stellar 
birth (Tc is the galaxy age and t(m) is the present stellar age), 
and f tp- agb {m) is the TP-AGB lifetime of a star with initial 
mass m. 

Let us now consider the typical mass range for the occur- 
rence of HBB, say from m = M HBB to m = M UB . The lower 

' J min U P 

limit corresponds to the minimum stellar mass necessary to get 
Tbce high enough to activate nuclear H-burning in the deep- 
est layers of the convective envelope, i.e. » lx 10 7 K for the 
p-p chains or « 2 x 10 7 K for the CNO cycle. Present stel- 
lar models (Romano et al. 2001) indicate that the production of 
lithium from AGB stars via the beryllium-transport mechanism 
requires M HBB * 3.O-3.5M as the metallicity increases in the 
range Z = 0.001 -0.02. 

The upper limit M up is the maximum initial mass for a star 
to develop an electron-degenerate C-O core, and its predicted 
value is sensitive to metallicity and model details, e.g. the treat- 
ment of convective boundaries. In most cases M up is expected 
to lie within 5 and 8 M . 

Relative to the mass range [M™ B ,M up ], what is the ex- 
pected fraction of AGB stars for which the present study could 
significantly alter the existing model predictions, i.e. in terms 
of evolutionary properties and chemical yields? 

To answer this question it is useful to calculate the cumula- 
tive distribution function 

rM 

J m hbb N(m) dm 

F(M) = P(m < M) = -£ (3) 

J m „bb N(m) dm 

where N(m) dm is given by Eq. ©. This latter gives the prob- 
ability that an AGB star has a progenitor with an initial mass 
m < M, with m varying within the range [M HBB , M up ]. 

For the sake of simplicity, a first estimate is obtained as- 
suming constant metallicity, constant SFR, ?tp-agb indepen- 
dent of the stellar mass, and adopting the Salpeter IMF {(p(m) oc 
m -235y jjj e resu it s are shown in Fig. [5] for three choices of 
M un and setting M HBB = 3.3 M . We see that in all cases 

up to nun u 

F(3.3M ) = and F(M up ) = 1 by construction. 

It turns out that F(4) = 0.33, 0.42, 0.53 for M up = 8,6,5M 
respectively, meaning that AGB stars with initial masses in 




3 4 5 6 7 8 

M/M Q 



Fig. 5. Cumulative distribution function defined by Eq. 
yielding the probability to have an AGB star with initial mass 
3.3 M Q < m < M. The distribution is normalised to the total 
number of AGB stars with masses in the range [M™ B ,M up ], 
where M up is set equal to 5, 6, and 8 M in the three plotted 
curves 

the interval [3.3, 4.0] M correspond to a sizable fraction of 
whole massive range. We remark that these are lower limits 
to the true values since larger fractions are expected when tak- 
ing the mass-dependence of ?tp-agb into account, as the TP- 
AGB lifetime should decrease with increasing stellar mass (e.g. 
Vassiliadis & Wood 1993). 

In summary, these simple calculations indicate that the im- 
pact of the variable molecular opacities should apply to an im- 
portant fraction of more massive AGB models. For our partic- 
ular class of models, with M up = 5 M , more than half of the 
stars experiencing HBB should be affected by the results of the 
present study. 

8. Cautionary remarks 

The results of this study depend on the assumed law for mass 
loss and its sensitiveness to changes in T e g . In our exploratory 
calculations we have adopted the formalism by Vassiliadis & 
Wood (1993) according to which, before the onset of the super- 
wind, the mass-loss rate scales exponentially with the pulsation 
period and P oc T^ M . This dependence on T e g derives from 
the underlying period-mass-radius for the fundamental mode 
(Wood 1990). We should remark that this latter relation is based 
on pulsation models for variable AGB stars with oxygen-rich 
envelope compositions. 

Due to the present lack of pulsation models suitably con- 
structed for C-stars, we have employed the same period-mass- 
radius relation for TP-AGB models regardless of the actual C/O 
ratio. As already mentioned, this fact implies that as soon as the 
third dredge-up makes C/O > 1, P is predicted to increase be- 
cause of the decrease in T e g caused by the enhanced molecular 
opacities. 

The physical soundness of this assumption could be veri- 
fied only with the aid of pulsation models for C-rich compo- 
sitions, using adequate /c var opacities. Preliminary calculations 
indicate that passing from C/O < 1 to C/0~ 2.5 - 3, the pul- 
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sation period actually almost doubles (Peter R. Wood, private 
communication). A detailed and quantitative investigation of 
this issue is left to a future analysis. 

Another source of uncertainty is related to the incapabil- 
ity of the our TP-AGB models, based on integrations of static 
envelopes, to account for the possible feedback between the 
third dredge-up and HBB as shown by Stancliffe et al. (2004). 
In brief, deep dredge-up episodes would actually increase 7b ce 
compared to models with shallower dredge-up, an effect which 
goes in the opposite direction to that due to the variable molec- 
ular opacities discussed in this work. In this context detailed 
AGB models with variable opacities are needed to fully inves- 
tigate the reciprocal feed-backs between the third dredge-up 
and HBB. 

9. Summary and conclusions 

This study has shown that hot-bottom burning may be weak- 
ened, quenched or even prevented as a consequence of the third 
dredge-up in intermediate-mass AGB stars. 

The effect of the third dredge-up against HBB is to be as- 
cribed to the response of the outer layers (i.e. convective en- 
velope and atmosphere) to the drastic change in the molecular 
opacities at the transition from C/0< 1 to C/0> 1 . Essentially 
two factors contribute, namely: 

- The sudden decrease of T e ff which accelerates the attain- 
ment of larger mass-loss rates, hence the reduction of the 
envelope mass; 

- The decrease of the temperature at the base of the convec- 
tive envelope, which weakens the nuclear reaction rates of 
the CNO-cycle. 

These circumstances are expected to take place only un- 
der specific conditions, namely i) in the early stages of the TP- 
AGB phase, the third dredge-up is able to convert the AGB star 
into a carbon star before hot-bottom burning becomes strong 
enough to efficiently destroy the newly dredged-up carbon in 
favour of nitrogen, and ii) the metallicity is high enough to al- 
low molecules to form in sufficient concentrations in the atmo- 
sphere. 

Our calculations indicate that such conditions are likely 
met in AGB models within a relatively small range of stellar 
masses, i.e. Ms3.0 - 4.OM and for metallicities Z ^ 0.001, 
while they should not apply to more massive models where any 
early transition to the C-rich domain is prevented by HBB it- 
self, and to more metal poor models for which molecules do 
not play an important role. 

Despite the narrowness of the involved mass range, it has 
been shown that the new results should affect a sizable fraction, 
say ^ 30 - 50%, of the whole group of massive AGB stars that 
could undergo HBB. 

The impact of such effect, so far neglected by full AGB 
models, may be large, so as to likely modify the current pre- 
dictions (e.g. nucleosynthesis, mass range, lifetimes, chemical 
yields) for the most massive carbon stars, that are predicted to 
populate the bright wing of the carbon star luminosity function. 
These luminous carbon stars are expected to exhibit the finger- 
prints of the HBB nucleosynthesis and they could explain the 



existence of particular groups of observed objects, such as the 
lithium-rich carbon stars (Smith et al. 1995), and possibly those 
among the J-type carbon stars that are found to be brighter than 
N-type stars (Morgan et al. 2003). Moreover, the new results 
could affect the predicted chemical abundances of planetary 
nebulae (PN) produced by intermediate-mass progenitors, i.e. 
by reducing the mass range of those that are usually associated 
to type I PNe (with measured over-abundances of He and N). 

Finally, our study indicates that the impact of variable opac- 
ities on the yields of the lowest metallicity massive AGB mod- 
els (Z < 0.001) is probably minor, due to the high dependence 
of molecular formation on temperature. 

Therefore, in order to perform a reliable comparison be- 
tween models and observations and to evaluate the role of these 
stars in the chemical enrichment of the interstellar medium, 
new AGB evolutionary calculations, including variable molec- 
ular opacities, are demanded. 
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